Symmetry breaking and multi-hump solitons in inhomogeneous gain landscapes 
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We address one-dimensional soliton formation in the cubic nonlinear medium with two-photon 
absorption and transversally inhomogeneous gain landscape consisting of a single or several am- 
plifying channels. Existence of the solitons requires certain threshold gain while the properties of 
solitons strongly depend on whether the number of the amplifying channels is odd or even. In the 
former case increase of the gain leads to a symmetry breaking, which occurs through the pitchfork 
bifurcation, and to emergence of a single or several co-existing stable asymmetric modes. In the 
case of even number of amplifying channels we have found only asymmetric stable states. 

PACS numbers: 



Emergence of localized nonlinear patterns supported 
by localized gain recently attracted increasing atten- 
tion. There have been reported stable one-dimcnsional 
(ID) structures in media described by the cubic complex 
Guinzburg-Landau equation with linear losses and either 
one [l| or two [H highly localized "hot spots". Stable 
solitons were also found in periodic lattices with active 
channels and nonlinear two-photon absorption Q, and 
at the surface of a periodic medium Q . Stable patterns 
may exist even in 2D settings, in layered structures of 
the planar waveguides ^ and in 2D lattices Q. 

Localized gain significantly changes the physics of 
emergent nonlinear patterns. In addition to standard 
constraints imposed by the balance between the dissipa- 
tion and gain, it introduces a new spatial scale - the 
width of the gain domain. This suggests a possibil- 
ity of existence of more sophisticated structures, than 
the simplest symmetric and/or anti-symmetric dissipa- 
tive solitons. In particular, it is of fundamental inter- 
est to exploit the phenomenon of the symmetry break- 
ing, which is known to be generic for conservative non- 
linear systems possessing a characteristic spatial scale. 
These are for example, systems governed by the nonlin- 
ear Schrodinger equation with a symmetric double-well 
potential tS, |9j . Self-trapping in one of the two identi- 
cal channels, has already been observed experimentally 
in nonlinear optics and in a Bose-Einstein conden- 
sate [lH. Symmetry breaking in a dual-core dissipative 
fiber with parametric gain and linear losses, which oc- 
curs due to interplay of the nonlinearity and coupling 
between two cores defined by the modulation of the lin- 
ear refractive index, was reported in 0] . Asymmetric dis- 
sipative solitons were reported only in 2D systems, e.g. 
in cubic-quintic Ginzburg-Landau equation |l2ll and in 
media with saturable gain and absorption [l3|. Here we 
report a principally different scenario of the symmetry 
breaking. It occurs in a medium without any conserva- 
tive potential or modulation of the refractive index, but 
having nonlinear dissipation and localized gain. 



Our setting is related to dissipative solitons, observed 
at the wavelength 1319 nm in self-focusing electri cally 
pumped waveguides fabricated on an InP substrate |14| . 
In such structures the two-photon absorption usually 
ranges from 10~^ to 10~^ cm/GW while localized gain 
can be implemented by using segmented strip-like elec- 
trodes or spatially-localized optical pumping. Solitons 
can be excited in 0.6 /im-thick planar guiding layers with 
a length of several millimeters, by input beams with typ- 
ical waists of a few micrometers, at gain levels about 70 
cm~^ [l4|. The considered model is also relevant for de- 



scription of Bosc-Einstcin condensates of quasiparticles 
in presence of nonresonant pumping [isj . 

We consider the propagation of laser radiation in a fo- 
cusing cubic medium with strong two-photon absorption 
and transversally inhomogeneous gain described by the 
equation for the dimensionless light field amplitude q: 
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+ ip^R{'n)q-\q\'^q-ia\q\'^q (1) 



where ^ and rj are the normalized longitudinal and trans- 
verse coordinates, respectively; pi > is the gain pa- 
rameter; R{r]) describes the transverse gain profile (its 
amplitude is normalized to one); a > is the strength 
of two-photon absorption. We consider gain landscapes 
containing an integer number n of periods of cos^ ij. For 
example, to model a gain with an odd number of ampli- 
fying channels we set R{'r]) = cos^ rj for jTyj < where 
rjn = mr/2, and R{ri) = for \ri\ > rjn, and vary pi, a 
and n. Assuming the characteristic transverse scale to be 
3 /im, we estimate the longitudinal scale (the diffraction 
length) to be ~ 170 //m at the wavelength of 1.32 /im. 
Then pi = 1 corresponds to the linear gain ^ 60 cm~ ^ 
and a = 1 corresponds to the two-photon absorption co- 
efficient « 0.017 cm/GW. The linear absorption is sup- 
posed to be compensated by the gain for \r]\ > rjn- 

Dissipative solitons of Eq. ([T]) can be searched in the 
form ij} = w^ifje^^^ , where h is the propagation constant, 
^(77) = Wr + iwi = u{rf)e^^^^^ is a complex amplitude. 
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with real, Wr, and imaginary, Wi, parts. The modulus u 
and phase 9 solve the equations 

bu=^+u'-£^^ j,, = 2p.,R{r,)u'-2au\ (2) 

where ^'(77) = can be referred to as a current density. 
We are interested in localized solutions with u,j — > 
at 1 77 1 00, which can be obtained for 6 > (more 
specifically with the exponentially decaying asymptotics 
M^e'^l"! andj ^e-4y2h|r,|)^ 

Localized modes of the system ([2]) form if a focusing 
nonlinearity counterbalances diffractive broadening, i.e. 
when the following equation 

bU=-^J uldi-i -J =^drj + J u'^di'i (3) 

where U = J u^drj is the energy flow, is satisfied and 
when the nonlinear losses integrally compensate the spa- 
tially inhomogeneous gain, i.e. when 

Pi j R{r])u^dT] — J u^dr] (4) 

The above formulas allow one to argue on possibility 
of existing of two different types of the modes at the 
same parameters of the system. Consider the limit of 
high amplification pi 00. Assuming that the solu- 
tion amplitude A ~ max{ii} grows and the width £ de- 
creases, the relation (j4|) suggests the scaling A ^ l/£ ^ 
^Jpija. This allows us to deduce from ^ the estimate 
\J « (ajpi) J u'^drj ^ y^pi/a, valid subject to the as- 
sumption that the soliton is maximum is placed exactly 
at 77 = where the pump has the maximal value, i.e. 
valid for a symmetric mode. For a crude guess of the pro- 
portionality coefficient in this estimate we use the ansatz 
u ~ i^/cosh(i'77) [that corresponds to neglecting the cur- 
rent j in the first of equations ([2]), what strictly speaking 
can be done only in the vicinity of 77 = 0]. This, allows us 
to obtain from for the symmetric mode ly = y/3pi/2a 
and thus U 2v = \J 6pi / a and b ^ v^/2 ~ ipi/Aa. 

However constraint (|4]) admits yet another scaling 
where a solution width grows with pi. Then for wide 
solutions, i.e. at € 770 = tt, from Eq. dH the relation 
Pi ~ aA?£ follows. On the other hand, now U ^ A^£, 
i.e. U ^ PijoL- Thus, unlike in the previous estimate, 
now we are restricted neither by the position of the max- 
imum of the mode, nor by the symmetry of its shape. 
Moreover, in the corresponding solution, the diffraction 
term u^^ ^ A/£'^ cannot be compensated by the Kerr 
nonlinearity v? ^ A? alone, and the role of the current 
distribution, i.e. of p /u^, becomes crucial (it reduces 
the impact of the Kerr nonlinearity). Notice that the 
major influence of the current occurs not at the origin 
(where for the symmetric solutions it is exactly zero) but 
at some intermediate point 77* defined by the the condi- 
tion jr](j]*) = 0. Thus if a solution with the suggested 



scaling exists, it should have asymmetric shape, with the 
maximum located in the vicinity of the point rj^ (at least 
in the limit pi — > 00). 

Further information about the maximal field ampli- 
tude A can be obtained from Eq. ([2]). Indeed, for rj > rjn 
the current is decaying, 7,, = — 2Q;7t^ < 0, and is di- 
rected outwards the gain domain: j > (since j tends 
to zero as 77 ^ 00). This means that maxima of \j\ are 
achieved at some points located inside the gain domain, 
i.e. |77*| < 77„. In such points the amplitude of the field 
is given by = piR{rj^,) / a (notice that 77* itself depends 
on the gain coefficient). Considering the symmetric one- 
hump mode in the case of one gain channel [i.e. when 
fimctions 71(77) ^'(77) feature only single maximum], 
one has the two maxima of j,, at ±77* . and hence j,, > 
in the interval |t;| < r;*. Therefore the amplitude of the 
field is bounded by the interval u1 < A? < piR/a. 

The above prediction of symmetric and asymmetric 
modes was confirmed in simulations [Fig. [T]. We ob- 
served that while the growth of zero background is sup- 
pressed at large rj, the light concentrates inside the am- 
plifying channels. Strictly speaking this feature is typi- 
cal for symmetric modes. A maximum of an asymmet- 
ric mode is shifted from the gain peak and the width of 
the mode grows with pi, according to the estimates pre- 
sented above. This broadening of the soliton leads to the 
situation where an appreciable part of the light energy 
concentrates outside the gain channel for large pi . 
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FIG. 1: (Color online) Symmetric (a) and asymmetric (b) one- 
hump solitons at n = 1, Pi — 3.5, a = 1.2, and (c) and (d) 
respectively profiles of symmetric and asymmetric two-liump 
solitons at 71 = 2, Pi = 2.5, a = 1.8. The modes in panels (a) 
and (b) correspond to circles in Figs.[3|a) and[3l^b). Hereafter 
all quantities are plotted in arbitrary dimensionless units. 
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The observed asymmetry in the field modulus remains 
relatively small for all considered pi (it is most clearly vis- 
ible in Wr^i distributions), but becomes more pronounced 
in systems with larger number of channels [see Fig. [2]. 
The number of humps of stable solitons coincides with 
the number of the amplifying channels. In Fig. [5] we show 
the symmetric and two co-existing asymmetric modes for 
the case of three amplifying channels (we obtained simi- 
lar solitons in landscapes with n up to 20). 
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FIG. 2: (Color online) Profiles of a symmetric three-hump 
soliton at n — 3, Pi = 2.5 [it corresponds to the circle in 
Fig.[3{c)] (a), and two coexisting asymmetric three- hump soli- 
tons at n — 3, Pi = 4.0, (b) and (c). (d) One- hump soliton 
at n — 3, Pi = 4.0. The dashed line schematically shows the 
gain landscape. In all cases a — 1.8. 

Asymmetric states in a system where gain landscape 
is symmetric and all other parameters are uniform, i.e. 
the symmetry breaking, is an unexpected result. Indeed, 
unlike in conservative systems, the understanding of the 
phenomenon cannot be related to the energetic argu- 
ments. Our system also does not allow for reduction 
to a simpler discrete model, as this happens, say, in the 
case of a double- well potential. Moreover, in our case the 
symmetry breaking occurs even for a single gain channel 
in contrast to conservative systems (where at least two 
potential minima are required). 

We performed numerical study of the whole branches 
of the solutions and studied their stability [Fig. [3]. In 
Fig. [3] (a), (b) for n = 1 we observe two branches of the 
solutions (notice, the propagation constant b is not a free 
parameter), one of them corresponding to the symmetric 
solitons, and another one, bifurcating form the symmet- 
ric branch at certain value pi = p^'', that corresponds 
to the asymmetric solutions (having smaller amplitudes 
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FIG. 3: (Color online) Energy flow (a) and propagation con- 
stant (b) vs Pi for symmetric "s" and asymmetric "a" one- 
hump solitons at n = 1, a = 1.2. Energy flow (c) and per- 
turbation growth rate (d) v& pi for symmetric three-hump 
solitons at n = 3, a = 1.8. The dashed lines indicate the 
borders of stability domains pi = p'y for symmetric modes. 



and larger widths as compared to the symmetric ones). 
The dependences U(j>i) and 6(pi) for the both branches 
well reproduce the estimates presented above. The linear 
stability analysis of the modes is performed by plugging 
in the perturbed field 9 = (w + we"^^) e*''* into Eq. ^ 
and performing linearization around w. For odd numbers 
of amplifying channels, exactly at the bifurcation point 
pf the branch of symmetric solutions looses its stabil- 
ity, while the stable asymmetric branch emerges [see the 
dashed lines in Figs. [3l[a) and (b)]. Since the asymmet- 
ric modes appear in pairs (corresponding to the left and 
right shifts of the maximum outwards the origin) at the 
point where the symmetric mode becomes unstable we 
deal with the pitchfork bifurcation. 

For small pi symmetric solitons broaden dramatically 
and may expand far beyond the region with gain (there 
is always a flow of energy outwards amplifying region). 
Increase of pi results in growth of the peak amplitude and 
progressive localization of the soliton inside the amplify- 
ing domains. According to the above estimates for n = 1 
and for sufficiently small a the energy flow and propaga- 
tion constant of a symmetric soliton arc monotonically 
increasing functions of pi [Figs.|3l(a) and (b)]. Note, that 
while for small values of a the symmetric one-hump soli- 
tons can be found even for pi — >■ 0, for moderate and 
high alpha values such solitons exist only above certain 
minimal value of gain coefficient [see Fig. d] (a)]. 

For even n the symmetric modes appear unstable in the 
whole domain of existence, and the only stable modes are 
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asymmetric ones. In this case the dependencies U{pi) for 
symmetric and asymmetric modes do not overlap and no 
bifurcations occur. Except for stabihty, other properties 
of modes supported by even and odd number of ampUfy- 
ing channels are similar. Gain with multiple amplifying 
channels also supports solitons with the number of humps 
smaller than the number of the channels [Fig. [2l^d)]. 

In dissipative multi-hump solitons both Wr and Wi 
change their signs in neighboring channels with gain, 
while the field amplitude u is nonzero even in the regions 
between the channels. This follows from Indeed, 
let us assume that at some point rj the field is zero, i.e. 
u(ji) = 0. Since uiji) is nonnegative, in the vicinity of fj 
we have: u(ry) = O [{t] — rff^ and = O ((ry — 77)*). Ex- 
panding u{rj) and j{7]) in the Taylor series in the vicinity 
of fj we find subsequently that all the expansion coeffi- 
cients are zero, what means that if u becomes zero at 
some point, u{ni) = and ^(jy) = 0. 

The critical gain, at which the bifurcation occurs, in- 
creases almost linearly with a, so that the domain of 
stability of symmetric solitons pf > Pi > expands 
with a [Fig. Sl^a)]. In the case of symmetric multi- 
hump solitons the energy flow grows with pi monoton- 
ically, except for the narrow region close to the thresh- 
old Pi = p'"™ below which no multi-hump solitons can 
be found [Fig. El^c)]. Multi-hump solitons arc stable in 
the region adjacent to p'"™, but increase of pi results in 
their destabilization at pi = p^^ [dashed line in Fig. [3jc)] . 
A typical dependence of perturbation growth rate dr on 
gain parameter for n = 3 is shown in Fig. El^d). The sta- 
bility domain of symmetric three-hump soliton expands 
almost linearly with increase of a [Fig.UJb)]. 
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FIG. 4: (Color online) The domain of existence {pi > p{°"') 
and stability domain [pi^ > Pi > Pi°™) for one- hump (a) 
and tliree-hump (b) solitons on the plane {a,pi). For small 
a < 0.6 symmetric one-hump solitons can be obtained even 
when Pi — ^ but for moderate and high nonlinear losses they 
exist only above certain minimal gain pi = p'°" 

Destabilization of symmetric multi-hump states is ac- 
companied by the appearance of several stable branches 
of asymmetric multi-hump solitons. Thus asymmetric 
modes depicted in Figs, ^h) and[2][c) that are both sta- 
ble and corresponding symmetric unstable three-hump 
mode (not shown) coexist for the same values of p^, a. 

With increase of the number of the gain channels the 



picture becomes even richer. When the number of chan- 
nels is odd, the number of asymmetric modes that can 
be stable all together for fixed Pi and a values increases. 
This feature indicates on the presence of several stable at- 
tractors (multistability) in multichannel landscapes. The 
critical value of the gain coeflicient (i.e. the bifurcation 
point) Pi^ also grows with n, reaching however certain 
saturation value. In particular, for a = 1.5 this value is 
about 2.59 and it is reached already at n = 7. 

Summarizing, we reported the symmetry breaking of 
dissipative soliton supported by a single or multiple am- 
plifying channels embedded in the cubic medium with 
nonlinear losses, which occurs through the loss of the sta- 
bility of the symmetric family at a point of the pitchfork 
bifurcation. Since inside the stability domains solitons 
are attractors with sufficiently large basin they can be 
excited with a variety of regular or noisy input patterns. 
For Gaussian inputs and single amplifying channel sta- 
tionary solitons may form already after propagation over 
20-30 diffraction lengths. In the case of multiple gain 
channels the bifurcation leads to appearance of several 
stable asymetric modes. While the bifurcation type is 
the same as one leading to appearance of the self-trapped 
states in a conservative double well potential , here we 
deal with pure dissipative phenomenon and the system 
does not possesses any characteristic scale related to its 
conservative part. 
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